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O ■ Abstract 

^ ! We present a combined measurement of -Rb = r(Z bb)/r(Z ^ hadrons) and 

I I the semileptonic branching ratio of b quarks in Z decays, Br(b iuX), using 

double-tag methods. Two analyses are performed on one million hadronic Z de- 
' cays collected in 1994 and 1995. The first analysis exploits the capabilities of the 

silicon microvertex detector. The tagging of b-events is based on the large impact 
^ ! parameter of tracks from weak b-decays with respect to the e^e~ collision point. 

^ I In the second analysis, a high-pf lepton tag is used to enhance the b-component in 

the sample and its momentum spectrum is used to constrain the model dependent 
uncertainties in the semileptonic b-decay. The analyses are combined in order to 
provide precise determinations of and Br(b iuX): 

i?b = 0.2174 ± 0.0015 (stat.) ± 0.0028 (sys.); 
Br(b ^ iuX) = (10.16 ± 0.13 (stat.) ± 0.30 (sys.))%. 



Submitted to European Physics Journal 



1 Introduction 



The Z partial width into b-quark pairs is a relevant parameter for precision studies of the 
Standard Model (SM) p. Due to the high mass of the top quark and the large top-bottom 
weak coupling, the process Z ^ bb receives sizable contributions from vertex corrections in the 
SM. The partial width is usually measured as its relative contribution to the Z hadronic width, 
i?b = r(Z^bb)/r(Z — ^ hadrons), since many experimental and theoretical uncertainties 
cancel when forming this ratio. A deviation from the predicted SM value 0.2158 ± 0.0002 
would point to the existence of additional vertex corrections and therefore would be a signal 
for new physics. 

The semileptonic branching ratio of b-hadrons, Br(b — > iuX), can be expressed as: 

Br(b - iuX) = = ^^^^ , (1) 

Tall 2rb^£,yX + Tb^Tj/X + Thad 

where Fb^^j^x is the semileptonic decay width, i being an electron or a muon, Fb^ri/X is the 
semileptonic component involving tau production and Fhad is the partial width from purely 
hadronic decays. Present measurements of Br(b iuX) at LEP [3] are in slight disagreement 
with model- independent measurements performed at the T{4S) jl]. Z decays into b-quarks 
allow for the presence of b-baryons, which have a lower semileptonic branching fraction [Hj. 
Measurements of the average number of charmed hadrons in b-decays, A^^c, provide an indirect 
estimate of the b ccs component. The measured values in the Nc — Br(b iuX) plane show 
that a discrepancy is still present between the T(45') and LEP results, as well as between the 
T(4S') results and theoretical predictions [H]. 

In this paper, we present a combined measurement of i?b and Br(b —>■ iuX) obtained with 
the L3 detector i7j at LEP. Features that distinguish the production of bb pairs from lighter 
quark production are: the long lifetime and hard fragmentation of b-flavoured hadrons, the 
large lepton momentum in semileptonic decays and the broad event shape caused by the large 
b-quark mass. The b-tagging method using lifetime information relies on data taken with the 
L3 Silicon Microvertex Detector (SMD) [Hj in 1994-1995. The method using the characteristic 
semileptonic decays of b-hadrons relies on the good lepton identification and lepton energy 
resolution of the L3 detector. It requires lepton candidates with high momentum along and 
transverse to the direction of the associated jet, caused by the hard fragmentation and high 
mass of the decaying hadron. 



2 Hadronic Event Selection 

Hadronic events are selected with criteria similar to the ones used for the measurement of the 
total hadronic cross section jHj. The basic requirements used there select 99.15% of all hadronic 
Z decays with a background of 0.15% from other processes. A few additional requirements 
ensure a good performance of the b-tagging techniques on the events. These are: 

• The number of reconstructed tracks must be larger than four. 
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• The polar angle of the thrust axis, 6't, reconstructed from calorimetric information, must 
be within a fiducial region delimited by the barrel part of the detector, | cos6'r| < 0.7. 

With these additional criteria, 970k hadronic events are selected from data, corresponding to 
an integrated luminosity of 71 pb^^. 

Hadronic Z decays are simulated using the JETSET (TU] generator and passed through a 
simulation of the L3 detector jTTj. The simulation takes into account the inefficiencies and reso- 
lutions of the different subdetectors as a function of time, weighted according to the integrated 
luminosity in data. 



3 Double- Tag Methods 



An event is first split into two hemispheres defined by the plane normal to the thrust axis. 
Separately for each hemisphere, a set of criteria is applied in order to significantly enhance 
the purity of b-events. A hemisphere satisfying those criteria is declared to be "tagged" . The 
number of tagged hemispheres, N^, is related to the total number of hadronic events, Nhad, by 
the following equation: 



2N> 



Rh^h + Rc^c + (l — -Rc — Rh) Cuds, (2) 



had 



where eb, Cc and euds are the tagging efficiencies for b, c and light quark hemispheres. The 
parameter i?b represents the ratio of cross sections (T(bb)/cr(hadrons). It is related to R^ by 
i?b = -Rb + 0.0003. The correction takes into account the contribution from photon exchange 
The shift is also present for i?c, but its effect on the measurement can be ignored for samples of 
reasonable b-purity. Therefore we will assume Rc = Rc in the following. The number of events 
with both hemispheres tagged, A^^, is given by: 

CbR^el + CcRcel + c^ds (l - i?c - ^b) eLs' (3) 



Nhad 

where the additional factors Cb, Cc and Cuds, called hemisphere correlation factors, quantify 
residual correlations between the two hemispheres, which lead to a deviation from the simple 
power law reduction of the efficiencies. 

Two parameters can be determined directly from data using the two experimentally mea- 
sured ratios: Nt/ {2Nhad) and Nu/Nhad', we choose to extract the parameters i?b and eb. The rel- 
ative rate of c production, Rc, is constrained to its experimental value -Rc = 0.1734±0.0048 jl2j . 
The efficiencies for quarks lighter than the b-quark, as well as the hemisphere correlation factors, 
are taken from the Monte Carlo (MC) simulation. However, only the factor Cb is relevant to 
the analysis for samples of good b purity (cc = Cuds = !)• Typical values of the MC parameters 
for high purity tags are ec ~ 2%, euds ~ 0.5% and Cb ~ 1. 
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If a different b-tagging algorithm is applied, three additional equations can be derived: 

Rh^h + Rc^'c + — Rc — Rh) Cuds) (4) 



Nt' 



§^ = c^R^e'^' + R.e',' + {l - R, - R^) e'J,, (5) 



2N, 



had 



Cb-RbebCb + Rc^c^c + (l — Rc — Rhj eudsCuds- (6) 



where the new efficiencies and correlations have similar meanings to the efficiencies and corre- 
lations defined for the tag t. The measurement can be performed by a global fit in which the 
values of the five ratios: Nt/{2 N^ad), Nu/Nhad, Nt>/{2 N^ad), Ntf/{2 N^ad): Nt^t'/Nhad are used 
to determine the values of R^, eb and e\^, while Cb, ec, ^uds, Cb, e^, e^jg, c'{, are constrained to the 
values obtained from the MC simulation with their statistical and systematic errors. 

The two tagging methods applied in this analysis are an impact parameter tag whose effi- 
ciency is denoted by e and a leptonic tag with efficiency e'. 



4 Impact Parameter Analysis 

4.1 Track and Primary Vertex Reconstruction 

The inner tracker of L3 reconstructs particle trajectories from hits in the two layers of double- 
sided silicon sensors of the SMD, up to 62 measurements in the central tracking chamber and 
two measurements in the Z chamber. These measurements are combined to obtain the five 
parameters characterising the trajectory, i.e. its curvature in the r — (p plane, its transverse 
distance of closest approach (DCA) to the vertex, its azimuthal angle at the DCA, its polar 
angle 6 and the Z coordinate at the DCA. Their covariance matrix is determined from the 
estimated single-point resolution function. The most important parameters for this analysis 
are the DCA and its error, aocA- Small biases in the DCA itself are removed by recalibrating 
the mean DCA value as a function of the azimuthal angle of the track and as a function of 
the track position inside a sector of the central tracking chamber The DCA error is 

recalibrated using tracks with a high probability of coming from the primary vertex. For these, 
the width of the DCA distribution for tracks with high momentum, where multiple scattering 
is negligible, is used to determine a factor that multiplies the calculated DCA error from the 
track fit. The factor is found to be close to one, in agreement with an analysis performed using 
high momentum tracks from e"'"e~ —>■ e"'"e~, fi~^fj^~ and t~^t~ events. Typical values for a^cA 
are 30 /im and 100 /xm for tracks with and without SMD information, respectively. In addition, 
the contribution from multiple scattering, not included in the error calculated in the track fit, 
is estimated from the dependence of the distribution width on transverse momentum. It is 
found that the additional multiple scattering error is 110/(p_LVsin^) fim for tracks with a hit 
in the inner SMD layer and 200/(p_LVsin6') for tracks without such hits, with p± measured 
in GeV. 

The average position of the LEP luminous region inside L3 is reconstructed using tracks 
collected in hadronic events. The position and its error are averaged over 200 consecutive 
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hadronic events, in order to follow drifts in the beam position. The result, called the beam- 
spot position, is used as a constraint in the reconstruction of the primary vertex in each event, 
weighted by the r.m.s. width of the luminous region in the horizontal (110 /im) and vertical 
(20 /im) directions. 

For the reconstruction of primary vertices, tracks are selected using the following criteria: 

• A track must consist of at least 20 hits in the central tracking chamber. 

• At least one hit in the inner layer of the SMD must be included in the track fit. 

• The DCA to the primary vertex has to be less than 1 mm. 

• The significance of the DCA, defined by the ratio of the DCA and its error, has to be less 
than five. 

• The transverse momentum of the track has to be greater than 150 MeV. 

The procedure uses an iterative method which starts from the beam-spot position as an initial 
estimate of the primary vertex position. At each step of the iteration, the vertex is calculated 
with all tracks selected for that step. If the probability of the vertex is less than 0.05, the 
track with the largest contribution to the is removed and the vertex is recalculated with 
the remaining tracks. This procedure is repeated until the probability of the vertex is at 
least 0.05 or only three tracks are left. At each step, the beam-spot position is used as a 
constraint. With this procedure, a primary vertex is reconstructed in 99.5% of the selected 
events. This efficiency is found to be independent of the quark flavour within one per mill, 
using MC simulation. 

The uncertainty on the vertex position depends on the azimuthal angle of the event thrust 
axis and on the number and quality of the tracks retained for its determination. The average 
uncertainty in the horizontal direction is 42 /xm for light-quark events and 77 /xm for b-events. 
The worse resolution for b-events is due to the unavoidable inclusion of b-decay tracks in the 
vertex determination. The vertical position of the primary vertex is known with high precision 
(20 fjLm) due to the small vertical width of the beam spot. 

4.2 Heavy Quark Tagging Using Impact Parameters 

The information from all tracks in a hemisphere is combined to form a discriminating variable, 
D, which describes the likelihood that all tracks come from the primary vertex. The sensitive 

single-track quantity used for constructing D is the impact parameter, defined as the absolute 
value of the track's DCA with a sign that is positive if the track intersects the direction of the 
associated jet in the direction of the jet's total momentum, negative if it intersects opposite to 
that direction. The angular resolution for the jet direction is 40 mrad. 

Tracks retained for the determination of D have to fulfil the following quality criteria: 



• The angle, 6j, between the track and its associated jet axis must satisfy cos6'j > 0.7. 
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• The track should have at least 30 hits spanning over a distance of at least 40 wires in the 
central tracking chamber. 



• The DCA to the primary vertex has to be less than 1.5 mm for tracks with SMD infor- 
mation. This cut is increased to 3 mm for tracks without SMD hits. 

• The angular separation of the track from the anode and cathode planes of the central 
tracking chamber, where the resolution is worse, must be more than 11 mrad. 

• If a track uses no hits from the SMD, at least 2 out of 8 hits from the inner portion of 
the central tracking chamber should be used in the track fit. 

The retained tracks are then grouped into different classes according to the pattern of the 
associated SMD hits. Each class corresponds to a different resolution function for the impact 
parameter measurement. The repartition among classes obtained in data is compared to the 
MC simulation. The proportions agree in absolute value to within a percent. 

The significance, s, is defined as the ratio of the impact parameter to its error. The total 
impact parameter error is composed of the error from the track fit, the multiple scattering 
contribution and the contribution of the primary vertex error, all determined according to the 
procedure described in Section H?!! 

The discriminant variable is constructed on the basis of a resolution function, R{x), which 
describes the probability that a track which comes from the primary vertex is measured to 
have an apparent impact parameter significance x. The probability, P, of finding a significance 
greater than the measured one, s, is given by: 

_ !:^Rix)dx 

The combined probability for the n tracks in the hemisphere is HILi P{si)- The probability, 
V{n), of measuring a value greater than OiLi is: 

i=l j=0 J- 

We define the discriminant variable as D = — \og^QV{n). A tagged hemisphere must have a 
value of D above some minimum value. This definition of the discriminant cut ensures that the 
amount of background from hemispheres which have tracks consistent with the primary vertex 
is 10~^, independent of n. 

The resolution function R is determined for each class by a fit to the significance distribution 
of all tracks with negative significance in data. The r.m.s. width of the significance distribution 
in each class is compatible with 1.0, but the distributions have substantial tails. Therefore, 
a model of the resolution function is constructed as a sum of two Gaussian functions and an 
exponential tail. The same resolution functions are used for data and MC. Special care is taken 
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to ensure that the simulation takes into account the multiple scattering and the time dependence 
of tracking chamber wire inefficiencies and wire resolutions as a function of the drift distance, 
SMD noise and SMD strip inefficiencies. Only small final adjustments are needed in order to 
reflect the behaviour of the data. The uncertainty on the adjustment will be used later for the 
determination of the systematic error due to resolution effects. Figure ^ shows the significance 
distribution for all tracks. For negative values of the significance, the good agreement between 
data and MC shows that the resolution effects are well understood. The positive part of the 
significance distribution is sensitive to the value of i?b- The data shows agreement with the 
MC distribution, which corresponds to a value of i?b = 0.217. 

The resolution function determined from data is used to calculate the track probability 
P{s). The distribution of the hemisphere discriminant D is shown in Figure |21 together with 
the MC expectation and its components, in terms of primary quark flavours. The agreement 
is satisfactory for the bulk, as well as the tail of the distribution, and the tagging power of the 
discriminant is clearly exhibited. 

The efficiency for b-tagging and the purity of the obtained sample can be varied by changing 
the cut on the discriminant variable. We obtain a b-tagging efficiency as a function of the sample 
purity as shown in Figure El for data and for the MC simulation. There is a residual difference 
between the efficiency observed in data with respect to the one predicted by MC, which never 
exceeds 2.5%. It is independent of the discriminant cut value over a wide range and it is 
consistent with the estimated statistical and systematic errors. For the discriminant cut used in 
the double-tag analysis {D > 2.3) we obtain ej^^*^ = (23.74±0.19 (stat.)±0.22 (sys.))%, whereas 
the MC estimate is e^*" = (24.21 ± 0.03 (stat.) ± 1.58 (sys.))%, where the systematic error is 
dominated by the b-physics modelling uncertainties. These uncertainties do not propagate to 
an error in i?b since the MC efficiency is not used in its determination. 

4.3 Systematic Errors 
4.3.1 Tracking Resolution 

The tracking resolution function is determined from data alone. Its statistical accuracy is 
such that it causes a negligible uncertainty on the measurement, but a wrong description in 
the MC simulation influences the values of auxiliary parameters like the efficiencies for lighter 
quarks and the hemisphere correlation factors. In order to estimate the systematic error due to 
uncertainties in tracking resolution, two MC samples are used, one with the final adjustment, 
MCfinah ciiid the other corresponding to the 1-sigma resolution uncertainty, MCi^-. The sample 
MCfinai is found to produce a stable value of i?b as a function of the discriminant cut within 
statistical and systematic errors. MCio- is defined by the change in the final adjustment that 
leads to a behaviour of i?b which is inconsistent with a constant value by one standard deviation 
of the observed fluctuations. The differences in ec, euds and Cb predicted by the two MC samples 
are propagated as an estimate of the error due to tracking resolution uncertainties. 
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4.3.2 Systematic Error from Background Modelling 

MC simulation is needed to determine euds, ^c, and Cb- For the charm efficiencies an accurate 
knowledge of production and decay properties of the charmed hadrons is important, since the 
different species, D°, D"*", Ds and Ac, have lifetimes varying in the range of 0.2 to 1.1 ps. 
Modelling uncertainties in euds arise from the residual contamination by light hadrons with 
long lifetime. Kg and A, as well as the probabilities for gluon splitting into bb and cc pairs. 
Modelling uncertainties of the b-hadron properties only influence the correlation factor Cb. They 
are estimated by varying the mean value of the B energy fraction < xe{^) >, the charged decay 
multiplicity and the average lifetime of b-hadrons. 

The variation of the parameters is performed following the suggestions of References [T^16j . 
The parameter ranges are listed in Tabled Tables |21 and El show the complete list of systematic 
uncertainties on ec and euds due to the propagation of these modelling uncertainties. The 
uncertainties on Cb are discussed in more detail in the following Section. 

4.3.3 Systematics from Hemisphere Correlations 

Systematic errors on Cb are due to uncertainties in the MC simulation. In addition to resolution 
and modelling effects, taken into account as described in Sections l4.3.1l and l4.3.2[ reconstruction 
algorithms and detector inhomogeneities may create correlations between the tagging efficiencies 
of both hemispheres. 

A possible source of correlation can be quantifled by choosing a variable A for each hemi- 
sphere which could be influenced by tagging the opposite hemisphere. For a particular cut on 
the hemisphere discriminant we then deflne three distributions: 

• The normalised distribution of A for all hemispheres, A^(A). 

• The single-hemisphere tagging efficiency as a function of A, e(A). 

• The normalised distribution of A in a co-tagged hemisphere, C(A). A co-tagged hemi- 
sphere is the one opposite to a tagged hemisphere, regardless of whether it is itself tagged. 

We then form a coefficient, c^, reflecting the correlation characterised by A for a particular 
discriminant cut: 

, ^ /6(A) C{X) dX 
Je{X)N{X)dX' 

The value Cb = 1 implies that there is no correlation between hemispheres from effects 
characterised by A. A value > 1 indicates a positive correlation, < 1 an anti-correlation. 

The analysis is performed on several candidate variables A. For a bb MC sample, the sum 
of the separate components is then compared to the total correlation factor Cb. A reasonable 
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Error source 


Variation 


Rc 


0.1734 ±0.0048 [12j 


Bottom fragmentation parameter: 




< XE{h) > 


0.709 ± 0.004 ITfJ 


Bottom decay parameters: 
B lifetimes 


1.55 ±0.05 ps 


B decay multiplicity 


4.955 ±0.062 [Ej 


Fractions in cc events: 




D+ 


0.233 ±0.027 [Ej 


Ds 


0.103 ±0.029 [IE! 


Ac 


0.063 ±0.028 [IE! 


Gluon splitting in cc events: 




g ^cc 
g ^hh 


(2.33 ±0.50)% [in! 
(0.269 ± 0.067)% [IH! 


Charm decay parameters: 




D" lifetime 


0.415 ± 0.004 ps [12J 


lifetime 


1.057 ±0.015 ps [12J 


Ds lifetime 


0.467 ±0.017 ps [12J 


Ac lifetime 


0.206 ±0.012 ps [12J 


D decay multiplicity: 
D*^ — s> prong 
D*^ — ^> 4 prong 
D*^ — s> 6 prong 


0.054 ±0.011 [IE! 
0.293 ±0.023 [IE! 
0.019 ±0.009 [IE! 


D"*" — s> 1 prong 


0.384 ± 0.023 [IE! 


5 prong 


0.075 ±0.015 [IE! 


Ds ^ 1 prong 


0.37 ±0.10 [IE! 


Ds ^ 5 prong 


0.21 ±0.11 [IE! 


D — >^ Kg multiplicity 


0.46 ± 0.06 TaHH! 


Charm fragmentation parameter: 




< xe{c) > 


0.484 ±0.008 


Fractions in uds events: 




KO and A 


JETSET ± 10% 


Gluon splitting in uds events: 




g ^cc 
g ^hh 


(2.33 ±0.50)% HE! 
(0.267 ±0.067)% fTH] 



Table 1: Variation of modelling parameters used for the determination of the systematic error 
in the impact parameter double-tag measurement. 
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agreement between them means that the relevant sources of correlation have been identified. 
The following sources are found to be relevant: 

• Angular effects: Inefficient regions of the detector can lead to correlations due to the back- 
to-back nature of hadronic events. This is estimated using A = cos 9 and A = 0, where 9 
and (j) are the polar and azimuthal angles of the most energetic jet in each hemisphere. 

• Vertex effects: Both hemispheres use the same primary vertex, which is determined as 
explained in Section 14.11 The main effect arises because the primary vertex fit may also 
include tracks coming from b-hadron decays (negative correlation). To quantify the effect, 
two independent primary vertices are constructed separately using the tracks assigned to 
each hemisphere. The variable A is taken to be the distance in the x — y plane between the 
vertex in each hemisphere and the overall event vertex. Its sign is given according to how 
far each hemisphere vertex moves when the beam-spot position constraint is removed: a 
positive sign is assigned to the hemisphere with the larger movement. 

• QCD effects: The presence of hard gluons in the event can influence the tagging efficiency 
of both hemispheres by taking energy away from the primary quarks (positive correlation) 
or, in an extreme case, by pushing both quarks into the same hemisphere (negative 
correlation). This effect is modifled by reconstruction and detector resolution effects. 
The signed event thrust, A = ±T, is used as a probe. A positive sign is assigned to the 
hemisphere with the higher energy jet. 

The dominant components are those due to gluon radiation and vertex bias. There is a 
remaining discrepancy between the linear sum of the correlation components due to the above 
three sources and the total observed correlation factor. That is expected, due to the interference 
between the sources considered and to additional sources of less relevance. For instance, and 
vertex effects are intrinsically correlated, since the primary vertex uncertainty is affected by 
the 0-dependent beam spot size. 

To test the quality of the MC correlation simulation, a 70% b-purity sample is selected in 
data by requiring the event discriminant to be greater than 1.5. The same cut is applied to 
the MC sample. The correlation coefficients are calculated for each of the three sources in data 
and MC. The differences are taken as the systematic uncertainties on the correlation term and 
propagated through to a systematic error on i?b- The complete list of systematic uncertainties 
on Cb is shown in Table EJ 

4.4 Results 

Systematic errors on ec, euds and Cb, are propagated into the double-tag i?b measurement. 
Figure E^i) shows the measured value of i?b as a function of the discriminant cut in the region 
around the minimum of the total uncertainty. The value is stable within the estimated statistical 
and systematic errors. Figure |3)d) shows the statistical and systematic errors on i?b in the same 
range. The minimum of the total error occurs at > 2.3, which defines the central value of 
our measurement. For this cut, 118817 tagged hemispheres and 11705 double-tagged events are 
selected from a sample of 968964 hadronic events. 



10 



Error source 


Acc 


MC Statistics 


0.02% 


Track Resolution 


0.02% 


D'*' fraction 


+0.10% 


Ds fraction 


+0.01% 


Ac fraction 


-0.03% 


D*^ lifetime 


+0.02% 


D"*" lifetime 


+0.01% 


Ds lifetime 


+0.01% 


Ac lifetime 


+0.01% 


Decay multiplicities: 




D''" 1-prong 


-0.03% 


D^ 5-prong 


0.00% 


D° 0-prong 


0.00% 


1)0 4-prong 


+0.02% 


I30 6-prong 


0.00% 


Ds 1-prong 


-0.02% 


Ds 5-prong 


-0.04% 


D — > Kg multiplicity 


+0.02% 


< xe{c) > 


+0.0D70 


g cc 


0.00% 


(7 — > bb 


0.00% 


Total 


0.14% 



Table 2: Error contributions to ec for a cut at > 2.3. The sign associated to a given error 
indicates the effect of a positive variation of the corresponding parameter in Table H 



Error source 


Aeuds 


MC Statistics 
Track Resolution 


0.00% 
0.01% 


K° and Hyperons 


+0.03% 


g ^ cc 
g^hh 


0.00% 
0.00% 


Total 


0.03% 



Table 3: Error contributions to euds for a cut at D > 2.3. The sign associated to a given error 
indicates the effect of a positive variation of the corresponding parameter in Table Ql 
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F/rroT RoiiTPP 




iviu statisxics 


U.UUoD 


Track Resolution 


n nnno 
U.UUUZ 


Vertex bias 


U.UUDo 




n finn? 


/7) POTFplfl'I'lOTlC: 
L-Wl 1 ClCLtlvJliiJ 


n nnriQ 


JTLdi U. g,lU.01i eililoolOli 


n nnriA 


B fragmentation 


+0.0014 


B lifetimes 


+0.0008 


B decay multiplicity 


+0.0004 


Total 


0.0067 



Table 4: Error contributions to Cb- The sign associated to a given error indicates the effect of 
a positive variation of the corresponding parameter in Table ^ 

The value of ec is estimated by MC to be ec = (3.05 ± 0.02 (stat.) ± 0.14 (sys.))%. A 
breakdown of the error is shown in Table |21 Since the individual charm lifetimes are measured 
very accurately, the fractions of the different species are the major error contributions. Among 
all the charmed hadrons the D+ properties lead to the dominant errors because it has the 
longest lifetime. 

The value of euds is estimated to be e^ds = (0.739 ± 0.004 (stat.) ± 0.035 (sys.))%. The 
different error contributions are shown in Table El The systematic error is dominated by the 
uncertainty in the rate of light-flavoured hadrons with long lifetime. 

The value of the correlation coefficient is Cb = 0.9717 + 0.0036 (stat.) + 0.0056 (sys.). The 
different error contributions are listed in Table |3] Primary vertex effects are the dominant 
source of uncertainty. 

The measured values of -Rb and are: 

i?b = 0.2173 + 0.0018 (stat.) + 0.0032 (sys.), (10) 
eb = (23.74 + 0.19 (stat.) +0.22 (sys.))%. (11) 

The detailed list of contributions to the error on i?b is given in Table The sources internal 
to L3 are separated from the ones in common with other experiments. 



5 Lepton Analysis 

The selection of leptons with high momentum and high transverse momentum with respect 
to the closest jet is applied to the set of events passing the hadronic selection. It requires 
a good performance of the central tracking detectors and a restriction to the angular region: 
I cos6't| < 0.7. 
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Ai?b from Internal Error Sources 


Mu statistics 


0.00092 


rtesoiULion 


u.uuuou 


A/p-rfp-v PTTPpft; on f^i. 


00125 


C7 eiiecLS un 


n nnnnfi 




n nnn9i 


Total Internal 


n nni p.p. 
U.UUioo 


Ai?b from External Error S 


ources 


Rc uncertainty 


-0.00094 


D"*" fraction 


-0.00128 


Dg fraction 


-0.00017 


Ac fraction 


+0.00045 


D"*" lifetime 


-0.00020 


D° lifetime 


-0.00013 


Ds lifetime 


-0.00012 


Ac lifetime 


-0.00008 


D"*" 1-prong decay multiplicity 


+0.00045 


D"*" 5-prong decay multiplicity 


-0.00006 


D*^ 0-prong decay multiplicity 


-0.00004 


D° 4-prong decay multiplicity 


-0.00024 


D° 6-prong decay multiplicity 


0.00000 


Ds 1-prong decay multiplicity 


+0.00028 


Ds 5-prong decay multiplicity 


+0.00055 


D — > Kg multiplicity 


-0.00025 


< xe{c) > 


-0.00086 


(? — cc in cc events 


-0.00001 


(7 — > bb in cc events 


-0.00003 


and A production 


-0.00176 


(yf — i> cc in uds events 


-0.00018 


g ho m uds events 


—0.00015 


B fragmentation < xeO^) > 


+0.00032 


B lifetimes 


+0.00018 


B decay multiplicity 


+0.00009 


Hard gluon fragmentation 


0.00010 


Total External 


0.00271 



Table 5: Internal and external error contributions to i?b for a cut aX D > 2.3. The errors 
from the D mesons fractions are propagated according to the correlation matrix defined in 
Reference jTHj. The sign associated to a given error indicates the effect of a positive variation 
of the corresponding parameter in Table Q] 
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5.1 Lepton Identification 

An electron is identified as a cluster in the electromagnetic calorimeter matched to a track in 
the central tracking chamber. The requirements are: 

• The cluster in the electromagnetic calorimeter must have more than 5 crystals with signal, 
more than 90% of the energy contained in a 3 x 3 matrix and a shower shape consistent 
with the deposit of an electromagnetic particle. 

• The cluster energy must be greater than 3 GeV. 

• The cluster must be matched in azimuthal angle (5 mrad) and in energy (within 4 standard 
deviations) to a track in the central tracking chamber. 

• The energy deposited in the hadron calorimeter in a cone of half opening angle of 7° 
around the cluster direction must be less than 3 GeV. 

A muon is identified as a track in the muon chambers which satisfies the following requirements: 

• The muon track is required to have track segments in at least two of the three layers 
of the chambers which measure the trajectory in the plane transverse to the beam. In 
addition, it must have at least one of the two hit segments in the muon chambers which 
measure longitudinal coordinates. 

• The measured muon momentum must be greater than 3 GeV. 

• The track must point back to the interaction vertex region within 3 standard deviations 
of the uncertainty in the extrapolation. This uncertainty takes into account multiple 
scattering and energy losses in the calorimeters. 

• The muon track must be matched to a track in the central tracking chamber. 

5.2 Heavy-Quark Tagging Using High-pt Leptons 

Hemispheres containing high-energy leptons are selected as b-candidates. Due to the hard 
b-fragmentation and the large b-mass, a higher b-purity sample is obtained as the transverse 
momentum of the lepton with respect to the b-jet, pt, increases. The distribution of this variable 
is shown in Figure El Very good agreement between data and simulation is observed. A pt cut 
of 1 GeV is used for the double-tag measurement. 

5.3 Systematic Errors 
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5.3.1 Lepton Selection Uncertainties 



Uncertainties in the lepton selection are expected if electromagnetic clusters, punch-through 
and tracking resolutions are not well reproduced in the simulation. Since the MC simulation 
takes into account the behaviour of the detector as a function of time, only a small uncertainty 
is expected. It is estimated by studying the statistical consistency of the measured R\, values 
as a function of the pt cut. 



5.3.2 Systematics from Background Modelling 



The charm efficiencies depend on the assumptions for the overall rate and properties of semilep- 
tonic c-decays. The uncertainty on these assumptions is taken into account by varying the 
semileptonic branching fraction, Br(c £i/X), the charm fragmentation, < xe{c) >, and the 
semileptonic decay model within the ranges suggested in Reference fHl- The parameter ranges 
are summarised in Table IHl The uncertainty on e^^g is estimated by selecting a sample enriched 
in light quarks. Both hemispheres in the event must satisfy a cut on the discriminant variable 
D < 0.5. In this sample a maximum discrepancy of 3% is found between data and MC for the 
number of leptons with momentum below 5 GeV. This discrepancy is translated into a relative 
error in e^^g. 



Error source 


Variation 


-Rc uncertainty 


0.1734 ±0.0048 ^12 


Charm fragmentation parameter: 




< Xe{c) > 


0.484 ±0.008 [in] 


Br(c ^ iuX) 


(9.8 ±0.5)% [16 


Semileptonic Decay Model c ^ i 


±ACCMM model 1 ™ 
-ACCMMmodel2 


Light quark uncertainties 


Data-MC comparisons 



Table 6: Variation of modelling parameters used for the determination of the systematic error 
in the lepton double-tag measurement. 



5.3.3 Systematics from Hemisphere Correlations 

The correlation between the semileptonic decay properties of the two hemispheres is expected 
to be very small. The value obtained for from the MC simulation is = 0.993 ± 0.009, 
consistent with no correlation. 

Evidence of non-negligible correlation effects have been looked for. The following quantity 
is defined: 

C = < "^"^ ^ , (12) 

<Xi> <X2> 
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where xi and X2 are values of physical variables in the two hemispheres of the event. We have 
used the lepton momentum and transverse momentum, since they are the relevant variables of 
the analysis. The values obtained for C are consistent with unity within the statistical errors, 
independent of the pt cut applied and of restrictions in the detector acceptance. Since no 
evidence of correlation is found, only the MC statistical uncertainty is assigned as a systematic 
error. 



5.4 Results 

For a Pt cut of 1 GeV, 49308 hemispheres are tagged and 1927 events are double-tagged on a 
sample of 968964 hadronic events. 

The value of is found to be = (1.70±0.01 (stat.)±0.11 (sys.))%. The error is dominated 
by the background modelling of semileptonic charm decays. The value of e^^g is determined to 
be e^ds = (0.362 ± 0.003 (stat.) ± 0.011 (sys.))%, and the value of the correlation coefficient is 
c{, = 0.993 ±0.009 (stat.). 

Taking into account statistical and systematic errors we obtain: 

i?b = 0.2116 ± 0.0050 (stat.) ± 0.0045 (sys.), (13) 
= (9.59 ±0.22 (stat.) ±0.11 (sys.))%. (14) 

The different contributions to the systematic error on i?b are shown in Table [7| 



Sources of systematics 


Error on i?b 


MC statistics 


0.0022 


Lepton Identification 


0.0005 


Rc 


-0.0013 


Br(c ^ iuX) 


-0.0028 


Semileptonic Decay Model 


-0.0022 


< xe{c) > 


-0.0006 


Light quark uncertainties 


-0.0007 


Total 


0.0045 



Table 7: Systematic error contributions to -Rb from the lepton double-tag measurement. The 
sign associated to a given error indicates the effect of a positive variation of the corresponding 
parameter in Tables ^ and IHl 



6 Measurement of i^b 

For an impact parameter discriminant cut of D > 2.3 and a lepton pt cut of 1 GeV, 9368 
events are selected with a lepton tag and an impact parameter tag in opposite hemispheres. 
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Ai?b from Internal Error Sources 


MC statistics 
Tracking Resolution 
Vertex effects on Cb 
6 effects on Cb 
effects on Cb 
Lepton identification 
Correlation between tags 


0.00081 
0.00042 
0.00079 
0.00003 
0.00013 
0.00010 
0.00079 


Total Internal 


0.00145 


Ai?b from External Error Sources 


Rc 


-0.00104 


D"*" fraction 
fraction 
Ac fraction 


-0.00100 
-0 00013 
+0.00035 


D"*" lifetime 
D*^ lifetime 
Dg lifetime 
Ac lifetime 


-0.00016 
-0.00010 
-0.00010 
-0.00006 


1-prong decay multiplicity 
5-prong decay multiplicity 
D*^ 0-prong decay multiplicity 
D° 4-prong decay multiplicity 
1)0 6-prong decay multiplicity 
Dg 1-prong decay multiplicity 
Dg 5-prong decay multiplicity 


+0.00035 
-0.00005 
-0.00003 
-0.00018 
0.00000 
+0.00022 
+0.00043 


D — i> Kg multiplicity 


-0.00019 


< xe{c) > 
Br(c IvX) 

Semileptonic Decay Model 


-0.00082 
-0.00068 
-0.00054 


(yf — i> cc in cc events 
(yf — s> bb in cc events 


-0.00001 
-0.00002 


and A in uds events 
Light quark uncertainties 


-0.00131 
-0.00007 


(yf — > cc in uds events 
(yf bb in uds events 


-0.00013 
-0.00011 


B fragmentation < ^^(b) > 
B lifetimes 
B decay multiplicity 
Hard gluon fragmentation 


+0.00020 
+0.00011 
+0.00006 
0.00007 


Total External 


0.00244 



Table 8: Internal and external error contributions to the i?b measurement. The errors from the 
D meson fractions are propagated according to the correlation matrix defined in Reference jl6j . 
Other errors in the table can be considered as uncorrelated. The sign associated to a given 
error indicates the effect of a positive variation of the corresponding parameter in Tables ^ and 

El 
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The correlation factor between the tags is = 1.004 ± 0.005 (stat.) ± 0.015 (sys.), where the 
systematic error takes into account detector and reconstruction effects. 

A combined fit using Equations 121101 for the impact parameter and lepton tags is performed. 
The MC input efficiencies and correlation factors with their uncertainties are included as con- 
straints in the fit. The final results are: 

i?b = 0.2174 ±0.0015 (stat.) ±0.0028 (sys.), (15) 
eb = (23.70 ±0.15 (stat.) ±0.19 (sys.))%, (16) 
e'^ = (9.34 ±0.07 (stat.) ±0.10 (sys.))%. (17) 

The statistical error of the measurement is determined by fixing the MC efficiencies and 
correlations to their fitted values. The systematic uncertainties are propagated according to 
the covariance matrix of the fit. Charm, light-quark and b-quark systematics are uncorrelated, 
except for the case of the charm fragmentation, which correlates the impact parameter and 
lepton tags. All systematic error contributions to the measurement are shown in Table |H1 

The measurement of -Rb is in agreement with the expectation from the Standard Model, 
j^SM ^ 0.2158 ± 0.0002 |2J, and is compatible with previous determinations |19j . 



7 Measurement of Br(b — > iiyX.) 

The efficiency e^, given by Equation (13 quantifies the fraction of high p,pt leptons in b-jets. 
Therefore, it is sensitive to the value of the semileptonic branching ratio of b-hadrons at LEP. 
In the reference MC with Br(b —>■ IvX) = 10.45% we find: 

g/^REF ^ ^g gg ^ Q_Q4 (^g^^^ ^) ^ Q_^^ (sys.))%. (18) 

The systematic error on is almost independent of b-quark model assumptions, since 
double-tag methods are used. The dependence on i?b is experimental, due to the simultaneous 
determination of both parameters (the statistical correlation is —0.72). The central value and 
the systematic uncertainty are determined for the set of parameters and variations shown in 
Table El The dominant systematic errors in are those due to detector inefficiencies and 

to semileptonic decay modelling uncertainties. They will be discussed in the next subsections. 

For different values of the branching ratio we find a linear dependence of e'^: 

e'^ = e'^^""^ + 0.5444 [Br(b ^ £uX) - Br(b ^ M)^^^]. (19) 

7.1 Detector Inefficiencies 

Time-dependent inefficiencies of the L3 detector are expected to be reproduced in the simulation 
within a few percent accuracy. In order to estimate possible extra sources of uncertainty not 
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taken into account in the MC, we have used large data samples of e"'"e~ e^e~ and e"'"e~ — 
fi~^fi~ . The samples are selected by requiring at least one lepton with high energy (> 0.35 ^/s) 
and applying the same cuts used for the sample described in Sectional except for the transverse 
momentum cut. The ratio of the number of one- lepton events to the number of two-lepton events 
is determined for both data and MC including backgrounds. The comparison suggests extra 
inefficiencies of (1.8 ± 0.1)% and (3.2 ± 0.1)% for the electron and muon samples, respectively. 
The errors are only statistical. These numbers are determined under the assumption that 
inefficiencies in both hemispheres of the event are uncorrelated. This is confirmed by the 
agreement at the percent level between the total cross sections measured in data and MC, 
which can be largely affected by simultaneous losses of the two leptons in an event. 

Additional checks at lower lepton energy (r ii^rT^e, e+e" e~^e~i~^i~) confirm the extra 
inefficiencies observed at high energy. Since the discrepancy is larger for muons, an additional 
cross check is performed. A sample is selected on the basis of the expected deposition in the 
hadron calorimeter for minimum ionising particles. The percentage of muons found in this 
sample is compared to MC, leading to an estimated extra inefficiency of (3.5 ± 0.5)%. When 
restricted to muons of low energy, we obtain (3.2 ± 1.2)%. 

As a consequence of this study the following relative losses have been added in the simula- 
tion: (1.8 ± 1.0)% for electrons and (3.2 ± 1.0)% for muons. The reference value in Equation 
UHl already takes into account these acceptance corrections. 



Parameter 


Variation 


Change in (%) 


e Detector efficiency 


±1% 


±0.037 


fj, Detector efficiency 


±1% 


±0.058 


Decay parameters in 6 — » £ 


Pf = 286 ± 35 MeV 


-0.080 


Decay scheme in 6 — £ 


±(ISGW** - ACCMM) 


-0.060 


< XE(h) > 


0.709 ±0.004 dH 


±0.032 


Br(b -> c -> £) 


(8.09 ±0.5)% 


-0.005 


Br(b ^ c ^ /) 


(1.66 ±0.40)% [lEI 


-0.010 


b-lifetime 


1.55 ±0.05 ps US] 


0.000 


Br(b ^T^i) 


(0.452 ± 0.074)% [in] 


-0.003 


Br(b ^ J ^ £) 


(0.07 ±0.02)% [ini 


-0.006 


Total 




0.11 



Table 9: Central values assumed for the MC determination of e[,. Systematic errors and ranges 
of variation are also shown. The sign associated to a given error indicates the effect of a positive 
variation of the parameter. 



7.2 Modelling Uncertainties 

The largest source of systematic error comes from the uncertainties in the modelling of the 
b — > iuX decay. The MC simulation uses a model which can be approximately described by 
the ACCMM [20] parametrisation with a c-quark mass rric = 1.67 GeV, a Fermi momentum 
inside the B meson Pf = 298 MeV, and a mass of the spectator quark rrisp = 150 MeV. In order 
to study the sensitivity to modelling assumptions, a sample enriched in b-quarks is selected from 
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data by applying an impact parameter discriminant cut of D > 1 on both event hemispheres 
and requiring the presence of at least one lepton with p > 3 GeV, pt > 1 GeV. The estimated 
b-purity of the 14929 selected events is 97.5%. 

Lepton hemispheres are classified in bins of (1 GeV x 1 GeV) in the {p,Pt) plane. The decay 
spectra in the context of the ACCMM model are generated as a function of rric and pf, and the 
MC {p,Pt) distributions are modified accordingly. A binned likelihood function C is defined: 

b' MC ^DATA 

= n" '^Bkri',' e-"."''--»). (20) 



where N^^"^^ and Nf^^{mc,Pf) are the number of events in the ith bin for data and MC, 
respectively. Only the shape is considered, that is, the total number of events in MC is nor- 
malised in order to agree with the total number in data. The data itself has no sensitivity 
for a simultaneous determination of rric and pf. On the other hand, different values of rric 
and Pf are found to reproduce the data with similar quality. Fixing the value of the c-quark 
mass to rric = 1-67 GeV leads to a minimum value oi pf = 286 ± 18 MeV, where the error is 
only statistical. Fixing the mass to rric = 1.60 GeV moves the central value to pj = 346 MeV, 
but leads to the same value of e^^^^. The quantity — 2 log(£/£Mc), where £mc is the value 
of the likelihood in MC at the minimum, should behave as a function. Its value is 1159.2 
for 1274 degrees of freedom. The previous fit is repeated without a cut on pt, with the result 
Pf = 273 ± 17 MeV and a of 1235.7 for 1340 degrees of freedom. A subsample with a harder 
lepton spectrum, p > 4 GeV,pi > 1 GeV, gives pf = 288 ± 20 MeV, with a similar fit qual- 
ity. If no impact parameter discriminant cut is applied, the value obtained is still statistically 
consistent: pf = 272 ± 13 MeV. 

Several systematic effects are studied. A 1% scale shift in the b-fragmentation parameter 
changes the value of pf by 14 MeV. A deterioration of the jet angular resolution by 20 mrad 
(half of the estimated resolution) leads to a shift of 20 MeV. Uncertainties in the momentum 
resolution give a smaller effect. We estimate a total systematic uncertainty of 30 MeV, yielding: 

Pf = 286 ± 18 (stat.) ± 30 (sys.) MeV. (21) 



Alternatively, the fit is performed in the context of the ISGW** model [21]. Using the D** 
fraction free parameter, we obtain: 

/(D**) = (24±4 (stat.) ±6 (sys.))%, (22) 



with a value -21og(£/£Mc) = 1150.3 for 1274 degrees of freedom. The ISGW** model with 
this central value leads to a slightly different value of e[,^^^. The difference is propagated as a 
systematic error due to the uncertainties in the modelling scheme. 

Finally, we perform a direct unfolding of the semileptonic decay momentum spectrum in 
the centre-of-mass frame of the b-hadron. For b-quarks, the original MC spectrum in the range 
— 2.4 GeV is divided into 12 bins. The relative contents of these bins are adjusted in order to 
make the simulated and the observed {p,Pt) spectra agree. The results are shown in Figure El 
There is good agreement with the optimal ACCMM, ISGW** spectra favoured by our data. 
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7.3 Results 



We use Equation El for the determination of Br(b — > iuX). Taking into account all the 
statistical and systematic errors on and Cb^^^, we obtain: 

Br(b ^ £uX) = (10.16 ± 0.13 (stat.) ± 0.30 (sys.))%. (23) 

Separate analyses, using only electron or muon tags, give consistent results within their 
statistical errors. The full breakdown of systematic errors is shown in Table The total 
systematic error due to modelling uncertainties on Br(b — > iuX.) is 0.2%. This error would 
increase to 0.3% if the /(D**) range of 11% to 32% is used jH]. However, this larger variation 
is excluded by our data at the 68% confidence level. The result represents an improvement 
on our previously published value and is consistent within errors, taking into account the 
different data samples, the different central values used in modelling and the fact that some 
systematic uncertainties are largely uncorrelated between the two measurements. 



1—' n T* n Try~i dl" d'v 

jr ai ameLei 


LxDiyu — ' tiyjvj \/o) 


MC statistics 


+0.07 


e detector efficiency 


-0.07 


fi detector efficiency 


-0.11 


Decay parameters in b ^ i 


+0.16 


Decay scheme in 6 — >■ £ 


+0.11 


< XE{h) > 


-0.06 


Br(b -^c^i) 


+0.01 


Br(b ^c^i) 


+0.02 


b-lifetime 


0.00 


Br(b ^T^i) 


+0.01 


Br(b ^ J ^ £) 


+0.01 


Uncertainty on Rc 


0.04 


Uncertainty on ec 


0.09 


Uncertainty on euds 


0.10 


Uncertainty in (c — > i) 


0.06 


Uncertainty in e'^^^ {i in uds) 


0.01 


Uncertainties on Cb, c'b, Cb 


0.08 


Total 


0.30 



Table 10: Systematic errors in the determination of Br(b — > iuX). The dominant uncertainties 
are the detector efficiency, the semileptonic decay model and the uncertainties coming from the 
use of the impact parameter tag in the determination of i?b- The sign associated to a given 
error indicates the effect of a positive variation of the corresponding parameter in Table El 

The semileptonic branching ratio for mesons containing a b-quark is extracted as follows: 

1 - /Ab 



21 



where f\^^ is the fraction of b-baryons in Z decays and Br(Ab — > iuX) is the semileptonic 
branching ratio in b-baryon decays. Assuming that all B-mesons have the same semileptonic 
branching ratio, our result is: 



Br(B^£z/X) = (10.47 ± 0.36 (stat.+sys.) (/aJ lals (Br(Ab £z/X)))% (25) 



where the last two errors are due to the uncertainties on the present measurements of /a^ = 
(10.ll|?)% [12] and Br(Ab ^ iuX) = (7.4 ± 1.1)% [5 . A similar conclusion is reached if the 
experimental values of the inclusive (rb) and b-baryon (ia^) lifetimes [T2j are used together 
with the assumption: Br(Ab £z/X)/rA^ = Br(b — ^z/X)/rb. 



8 Conclusion 

Using double-tag methods on one million hadronic Z events collected in 1994 and 1995, we 
determine the values of -Rb = r(Z —* bb)/r(Z — hadrons) and of the semileptonic branching 
ratio of b-quarks, Br(b iuX). Two methods of tagging the presence of b-quarks are used. 
The first method exploits the capabilities of the silicon microvertex detector and is based on the 
large impact parameter of tracks from weak b-decays with respect to the e"'"e~ collision point. 
In the second method, a high-p^ lepton tag is used to enrich the b-content of the sample and 
its {p,Pt) spectrum is used to constrain the model dependent uncertainties in the semileptonic 
b-decay. We measure: 



The measurement of i?b agrees with the Standard Model expectations and with previous 
determinations The measurement of Br(b — > £z/X) is in agreement with measurements of 
the semileptonic branching ratio of B mesons at the T[4S) 0. 
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Figure 1: Distribution of the significance. Tfie points are tfie data and tfie fiistogram is tfie 
MC prediction for a value of i?b = 0.217. 
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Figure 2: Distribution of the discriminant variable, D, in data compared with the MC pre- 
diction. The flavour composition of the hadronic sample is also shown. Good agreement is 
observed. The higher b-tagging power for large values of the discriminant is exhibited. 
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• Purity 

o MC efficiency 

o Data efficiency 
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Figure 3: Efficiency and purity of the sample obtained as a function of the cut on the discrim- 
inant variable, D. For a discriminant cut a,t D > 2.3, the purity is 0.843. 
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Figure 4: a) Value of i?b obtained as a function of the cut on the impact parameter discriminant, 
D. b) Statistical, systematic and total errors for i?b as a function of the discriminant cut. The 
uncertainty on Rc, ARc = 0.0048, is taken from Reference jl2j . 
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Figure 5: Distributions of the transverse momentum with respect to the closest jet for electrons 
a) and muons b). The transverse momentum of the lepton ^ is denoted by p^. The histogram 
shows the flavour composition of the MC sample. 
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• Unfolded spectrum 
— ACCMM, Pf=286 MeV 
■■■ ISGW, f(D**)=24 % 
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Figure 6: Spectrum of the lepton momentum, pi in the centre-of-mass frame of the semileptonic 
decaying b-hadron. The points and error bars are obtained by unfolding of the {p, pt) spectrum 
observed in data. All points are statistically correlated. The histograms correspond to the 
optimal spectra favoured by the data in the context of ACCMM and ISGW** models. Overall 
consistency is observed. 



32 



